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It is known that classical electromagnetic radiation at a frequency in resonance with 
energy splittings of atoms in a dielectric medium can be described using the classical sine- 
Gordon theory. In this paper we quantize the electromagnetic field and compute quantum 
corrections to the classical results by using known results from the sine-Gordon quantum 
field theory. 
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1. Introduction 



The importance of integrable non-linear partial differential equations in classical non- 
linear optics has been recognized for some time. The most well-known example concerns 
weakly non-linear dielectric media, where only the first non-linear susceptibility is consid- 
ered important. In this situation, the envelope of the electric field satisfies the non-linear 
Schrodinger equation. This was first understood theoretically by Hasegawa and Tappert[|l|]. 
The solitons predicted in [|l| were observed experimentally in 0. The occurence of these 
classical solitons in common optical fibers promises to revolutionize high-speed telecom- 
munications. 

Our interest in this subject concerns the possibly interesting quantum effects which 
arise when the electromagnetic field is quantized. The latter quantization amounts to 
studying an interacting quantum field theory in the classical variables which satisfy the 
non-linear differential equation. These quantum integrable models have been studied ex- 
tensively, and many of their properties have already been exactly computed. A priori, one 
expects quantum effects to be small for macroscopically large objects such as solitonJ^. 
Nevertheless, using known exact results from the quantum non-linear Schrodinger theory, 
quantum effects have been predicted and measured 

It is well-known that the non-linear dielectric susceptibilities are enhanced when the 
radiation is in resonance with the energy splitting of quantum states of the atoms of 
the sample. Near resonance, one is no longer in the weakly non-linear regime (higher 
susceptibilities involving higher powers of the electric field become as important) and the 
physics is no longer well described by the non-linear Schrodinger equation. Remarkably, 
as was understood and demonstrated experimentally by McCall and Hahn[0, the system 
is well described classically by another famous integrable equation, the sine-Gordon (SG) 
equation. 

In this paper, we study the quantum effects which arise when electric fields in res- 
onance with a dielectric medium are quantized by using known exact results for the SG 
quantum field theory. Though these quantum effects are theoretically interesting, as we 
will show, they are unfortunately probably too small to be measurable at the present time. 

^ In fiber optic systems the soliton consists of a cluster of 10^ or more photons. 
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2. Classical Theory 

In this section we review the manner in which the classical sine-Gordon equation arises 
in resonant dielectric media 0. 

We consider electromagnetic radiation of frequency u propagating through a collection 
of atoms, where the frequency uj is in resonance with an energy splitting of the atomic 
states. For simplicity, we suppose each atom is a two state system described by the 
hamiltonian Hq with the following eigenstates: Hq\iPi) = —^huo\tpi), -ffol'02) = |^'^o|'*/'2), 
such that hujQ is the energy difference of the two states. In the presence of radiation, the 
atomic hamiltonian is 

Hatom=Ho-p-E, (2.1) 

where p = e is the electric dipole moment operator. 

We assume the radiation is propagating in the z direction, and E = nE{z,t), where 
n - ^ = 0. The only non-zero matrix elements of the operator p ■ E can be parameterized 
as follows: 

{i;2\p-E\i;i)=pE{z,t)e-'", (2.2) 

where p and a are constants which depend on the atom in question. (We have assumed 
spherical symmetry.) It is convenient to introduce the Pauli matrices Ui, and write the 
hamiltonian as 

Hatom = —hujoas - E{z,t){piai +P2(r2), (2.3) 

where pi — pcosct, p2 = psina. 

The dynamics of the system is determined by Maxwell's equations, 

(d! - ^df^ E{z,t) = ^ dfP{z,t), (2.4) 

where = c^/eo, eo is the ambient dielectric constant, and P = nP{z^t) is the dipole 
moment per unit volume. The latter polarization can be expressed in terms of the expec- 
tations of the Pauli spin matrices (cXi) = ('i/'|a"i l?/^), where IV') is the atomic wavefunction. 
Namely, 

P = nn{pi{ai)+p2{a2)), (2.5) 
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where n is the number of atoms per unit volumeB. Thus, in addition to the MaxweU equa- 



tion (|2.4|) , one has dynamical equations for the polarization P[z,t) which are determined 
by Schrodinger's equation for the atom: 

ihdt^cJi) = {ip\[ai,Hatom]\'>P)- (2.6) 
The latter can be expressed as 

dM) = Y.e,jkVj {ak), (2.7) 

where e is the completely antisymmetric tensor with £123 = 1, and 

2E{z,t) 2E{z,t) 

Vi = pi, V2 = P2, V3 = ujQ. (2.8) 

n n 

To summarize, the dynamics is determined from the coupled equations of motion 

( p.4| ) and ( |2.7| ), wherein the atoms are treated quantum mechanically and the radiation is 

classical. 

Let E{z,t) = S{z, t) cos{u!t — kz) where uj/k = c and £{z,t) is the envelope of the 
electric field. We will assume the envelope is slowly varying in comparison to the harmonic 
oscillations: dt£ « uj£, d^S << kS. In this approximation, one finds 



d, + ^dt]S{z,t) 



sin{ut -kz). (2.9) 



Let us define (cr||), {a±) as follows 

{an) = (cTi) cos(a;t — kz + a) + ((72) sin(a;t — kz + a) 

(2.10) 

{a±) = — (cTi) sin(u;t — kz + a) + ((T2) cos(c<;t — kz + a). 

We make the further approximation that cos2(a;t — kz + a) terms in the equations of 
motion for dt{(Ji) can be dropped in comparison to cos(c<;t — kz + a) (and similarly for the 
sine terms)i. (This amounts to replacing sin^ [ut — kz + a), cos^ {ut — kz + a) by 1/2) . One 
then finds 

dt{a\\) = {u - ioo) {cri_) (2.11a) 
dt{ai_) = (^3) + {coo - u;){a\\) (2.116) 

dt{as) = ^^p{a^). (2.11c) 



^ To be more precise, {ai) here represents average over many atoms in a small volume, and is 
thus a continuous field depending on z,t. 

^ It can be shown in perturbation theory that this is a good approximation at or near resonance. 
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Finally, equations ( p.4|) and ( |2.9| ), upon making an approximation analagous to the slowly 
varying envelope on the RHS of lead to 



{d, + ldt)S{z,t) 
c 



27r 

sin(a;t — kz) = —cnup (sin(a;t — kz){a±) — cos{ujt — kz){a\\)) . 

(2.12) 

In order to solve these equations, note that dt{J2i{^i) i'^i)) — 0- Thus, if the atoms 
start out in their ground state with (era) = 1/2, one has J^ii'^i)'^ — 1/4 for all time. 
On resonance, when uj = ujq, eq. (|2.11a| ) implies (a||) = for all time. The constraint 



^j(ai)^ = 1/4 can be imposed with the following parameterization: 

{ai_) = ^ sm{f3,i<p{z, t)), (as) = ^ cos(/3c/(/)(z, t)). (2.13) 

The parameter Pel is arbitrary at this stage, but will be fixed in the next section. Equations 
( ^.116, cD now imply 



dt<P = --fT ^{z,t). (2.14) 
Inserting ( p.l4| ) into ( |2.12|) and defining x = 2z — ct, one obtains the SG equation: 

{d^-c'dl)<f>=-^smiPei<P), (2.15) 



where 



2 ^Tinp'^u 



3. Quantum Effects 

In this section we proceed to quantize the electromagnetic field. In order to do this, 
the SG field (p must be properly normalized such that the energy of soliton solutions 
corresponds to the physical energy; this amounts to properly fixing the constant f3ci- 

The action which gives the classical SG equation of motion is 

^SG = \j dxdt Q(c>i(^)2 - y (c>,(/>)2 + ^ cos(/3cZ(/>)) . (3.1) 
On the other hand, the properly normalized Maxwell action is 

^Maxwell = d^xdt (^-If^^F'^^ + =^ j d^xdt ^{dtAf + .... (3.2) 



where the vector potential is A = nA and as usual E = —^dfA. (One can show that Aq can 
be set to zero.) To normalize the field 4>, one needs only compare the kinetic terms in ( |3.1j ) 
and (^^)- The dimensional reduction is made by assuming simply that A is independent 
of y, z and / dydz = A, where A is an efi'ective cross-sectional area perpendicular to the 
direction of propagation (as in the cross-sectional area of a fiber). From ( |2.14|) one has 

5'Maxwell = ^ (~~) j ^^^^ + •••• (3-3) 

Comparing with (pT|), one fixes = 2p'^ / {Ah^c) . 

Finally, as is conventionally done, we rescale (p — > V^^, so that Ssg/^ takes the form 
(pT|)with (3ci replaced by /3, where 

The constant (3 is the conventional dimensionless coupling constant in the quantum SG 
theory, and is allowed to be in the range < < Stt, where 0^ = 8n corresponds to a 
phase transition[§. The limit where the radiation is classical but the atom is still quantum 
mechanical corresponds to the limit /3 ^ 0. 

The magnitude of the quantum corrections to classical results is determined by the 
parameter (3'^ /Sir. One can obtain an order of magnitude estimate by taking p ^ ei^bohr, 
where -Rbohr is the Bohr radius. One finds 

f^lO--^, (3.5) 
87r A 

where A is in cm^. Disappointingly, for realistic A, the constant (3'^/8n is exceedingly 
small. For the remainder of this paper we describe quantum effects that are at least in 
principle measureable, though in practice probably too small to observe. 

The classical soliton solutions to the SG equation are characterized by a topological 
charge T = ±1, where 

T=^((/)(a: = oo)-(^(x = -oo)). 

Solitons of either charge correspond to solutions where at fixed z the atoms in the far past 
are in their ground state, and are all in their excited state at some intermediate time: 
(o"3)t=±oo = 1/2. These classical solitons have been observed experimentally in 0]. What 



distinguishes solitons (T = 1) from antisolitons (T = —1) is the sign of the envelope of the 
electric field. From the known classical soliton solutions and ( |2.14| ) one finds for T = ±1, 



^(x,t) = ±^j|±^ cosh 4^11 . (3.6) 



P y C-V y yy/^ _y2 

Thus the electric fields for the soliton versus the antisoliton are out of phase by tt. 

The classical soliton mass is given by Ms = 8m//3^, where m = ^fi/l?. From 
one can express /3^/87r in terms of the classical soliton mass: 

Thus, 



Stt A^^ \-t^atom 

where A^-y = fiuj/MgC^ roughly corresponds classically to the number of photons that 
comprise the soliton, Ac is the Compton wavelength of the soliton, and Latom = ^/nA 
is the inter-atomic spacing. The above equation summarizes where one expects quantum 
effects to be important: when the soliton is comprised of small numbers of photons, or 
when the Compton wavelength is large compared to the space between the atoms. 

We now derive quantum corrections to the frequency and density dependence of the 
soliton mass. Classically, Ms oc ^Jum. In the quantum theory, short distance singularities 
are removed by suitably normal ordering the cos(/3(/)) potential^]: 

-r^ cos(/3(/)) A : cos(/3(/)) : 

P 

The anomalous scaling dimension of the operator : cos(/3(/)) : is /3^/47r, so that A has mass 

871- ' 

M, = c(/3)(v^)^/^^-i^) (3.9) 



dimension 2(1 — ^). Therefore, 



for some constant c(/9). The latter constant was computed exactly in [Q. Since A oc un 
one finds 



oc v^(^l + ^log(wn) +0(/3*)J . (3.10) 

The classical scattering matrix for the solitons has been computed by comparing N- 
soliton solutions in the far past and far future[0. The exact quantum S-matrix is also 
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known 1 11]. Since (3 /9>ti is small here, the quantum corrections to classical scattering are 



most easily determined by incorporating the one-loop corrections to the classical scattering, 
which amounts to the replacement = /3^/(l — |^). The necessary formulas can be 

found in the above papers. 

The most interesting aspect of the quantum scattering of solitons is that the so- 
called reflection amplitude is a purely quantum effect, analagous to barrier penetration 
|l'4 j. Namely, one considers an in-state consisting of a soliton of momentum pi and an 
antisoliton of momentum p2 which scatters into an out state where the momenta pi and p2 
are interchanged. The T = -fl soliton is thus reflected back with momentum equal to that 
of the incoming T = — 1 antisoliton. In the semi-classical approximation, this reflection 
amplitude is 

where 

^ (3.12) 



S {6) — exp ( — / drj log 

,7 Jo 



and 9 = 9i — 62, where pi^2 = Mg sinh^i^2- Since I/7 ~ l//3^ is very large, the oscillatory 
factor exp(167r^i/7) — 1 in ( |3.11D will make the detection of reflection processes difficult, 
since even in a small range of this factor averages to zero. 
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